Interchain versus intrachain hole transmission through desoxyribonucleic acid 

molecular wires 
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We present a methodology for computing the current-voltage response of a molecular wire within 
the Landauer-Buttiker formalism based upon transforming the cummulative transmission probabil- 
ity into an eigenvalue problem. The method is extremely simple to apply since does not involve 
construction of the molecular Greens function, and hence avoids the use of complex integration con- 
tours to avoid poles. We use this method to study the effect of base-pair sequence on the conductivity 
of holes in DNA chains containing A-T bridges between guanine chains. Our results indicate that se- 
quence plays a substantial role in ballistic transport via tunneling resonances tuned by sequence and 
interchain interactions. We also find that ballistic transport is dominated by intrachain transport 
and that hole transmission is insensitive to interchain fluctuations. 



I. INTRODUCTION 

Advances in synthetic technique and nano-scale device 
fabrication have facilitated the creation of prototypical 
electronic logic components in which the length scale of 
the device is on the order of the de Broglie wavelength 
of an electron. The prospect of using such "molecular 
electronic" components in technological applications is 
a strong driving force in the emerging science and en- 
gineering of nanotechnology. The real promise of nan- 
otechnology relies upon on our ability to predictably ma- 
nipulate matter on the molecular length scale. In partic- 
ular, for molecular electronic materials, the goal is to be 
able to fine tune the quantum mechanical energy states 
so that the properties of the individual quantum states 
is expressed whatever macroscopic device we are trying 
to engineering. The challenge here is then to build the 
interconnects between the individual molecules and the 
"outside world" so that one can reliably address and ac- 
cess the molecular-scale components. 

The sequence dependence of charge transport in DNA 
is a topic of considerable interest from both experimen- 
tal and theoretical viewpoints. From a design stand- 
point, DNA has a number of desirable properties. It 
possess a high degree of of site specific binding between 
single strands of DNA and its related self-assembly prop- 
erties. Furthermore, one can synthesize essentially any 
sequence desired and hence potentially tune its proper- 
ties in a controllable way. Such characteristics have made 
DNA an ideal candidate for incorporation into molec- 
ular scale electronic devices. With such motivation, a 
number of studies of the charge-transport properties of 
DNA have been performed 0HHIlllli,|10 and 
there is an emerging consensus regarding the mechanism 
for single-electron transport in DNA chains, [fj However, 
the DNA's intrinsic conductivity remains highly contro- 
versial. Depending upon the experiment. DNA is a.) an 
insulator at room temperture |10l llll 1 1 2L 1 1 3| , lx) a wide- 
band gap semiconductor at all temperatures 0,0], c.) 
Ohmic at room temperature and an insulator at low tem- 



pertures [H 111 IrlJTl llH I20L EH . or d.) metallic down 
to low tempertures [22| with induced superconductivity. 
A recent review by Endres, Cox, and Singh nicely sum- 
marizes these varied experiments. |23| 

A number of studies have been performed at various 
levels of theory with the aim of computing the current- 
voltage (IV) response of a molecular wire. 

HE H^, m, US, m, B EE HE m, l39l The studies 

differ widely in terms of the level of treatment of the 
electronic structure, the coupling between the molecule 
and the semi-infinite leads, and the change of electro- 
static potential due to bias across the system. Both first 
principles and semi-empirical treatments have been used, 
with the former being restricted to smaller scale systems. 

We consider the case where the terminal G's and G — C 
pairs are anchored to opposing semi-infinite continua cor- 
responding to the metal contacts. Similar models have 
been considered by Jortner and Bixon [4|, |8( consider- 
ing the canonical transfer rate from to G/ sepa- 
rated by (T — A) n bridges. Here, we consider the mi- 
crocanonical transmission probability at a given scat- 
tering energy, T(E). From this we can compute the 
transmission current as a function of applied voltage bias 

via|ll in m m m m EE S3 



W = | / T(E,V)(f L (E- 



V/2) + f R (E + V/2))dE(l) 



where fh,R are the Fermi distribution functions for the 
left and right hand electrodes and T(E, V) is the trans- 
mission probability at a given energy E and applied 
voltage bias V. We examine how permuting base-pairs 
changes the transmission probabilities and hence the 
current-voltage response. Aside from differences in pa- 
rameterization of the models, the results presented herein 
parallel those of Roche et al. who considered the the 
scaling behavior of the intrinsic conductivity of DNA 
due to length and sequence correlations. [42, Hlj Here 
we compare intra- vs. interchain hole transport whereas 
the work by Roche et al considered considered the DNA 
strand as a single quasi-one dimensional strand. 
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5' — G — T — T — G — 3' 

A) \ X 

3' — C — 'A' — A' — C — 5' 



5' — G — A — A — G — 3' 

B) 

3' — C — T — T — C — 5' 



5' — G — A — T — G — 3' 

C) X X 

3' — C — T — 'A' — C — 5' 



5' — G — T — A — G — 3' 
D) \, '-;<• 

3' — C — 'A' — 'T — C — 5' 

FIG. 1: Intrachain and interchain coupling pathways for hole 
transport in DNA sequences considered. Solid connecting 
lines correspond to intrachain hopping matrix-elements and 
dashed lined correspond to non-zero interchain hopping ma- 
trix elements. 



Here present a simple approach for embedding a 
generic quantum chemical problem into a scattering cal- 
culation for the purposes of computing the current- 
voltage response of molecular wire. Our approach uses 
a idea originally developed by Manthe and Miller for 
computing the canonical thermal rate constant for bi- 
molecular reactions in reactive scattering. The orig- 
inal idea was to apply absorbing boundary conditions 
about the transition state and use these to artificially 
impose the correct boundary conditions on the Greens 
functions used to construct the transmission matrix di- 
rectly, thus, avoiding the explicit construction of the scat- 
tering matrix. In the case at hand, we impose absorbing 
boundary conditions through the use of a complex self- 
energy term for atoms in our system in direct contact 
with semi-infinite leads. 



II. HOLE TRANSMISSION THROUGH MODEL 
DNA SEQUENCES 

A. Tight-binding model 

In what follows, we consider the bridging molecule to 
be a single DNA strand with matched nucleotide pairs 
attached to the left and right metal clusters at the ter- 
minal sites as shown schematically in Fig. [21 Guanine, G, 
has the lowest oxidation potential of the four nucleotides 
and we consider the case where an injected hole travels 
along the DNA chain hopping on the sites containing G. 



As the hole propagates along the chain, sites containing 
A-T pairs act as potential barrier between guanine sites 
that must be overcome by quantum mechanical tunnel- 
ing. Between guanine sites, hole transfer is taken to be 
effectively barrierless with no net change in free energy as 
the hole hops from G n to G n ±\. Since scattering dynam- 
ics is dominated entirely by what happens when the par- 
ticle encounters a potential change, the terminal G units 
can be considered to be effective sites linking a short seg- 
ment containing (T—A) bridges to the asymptotic region. 
When there is a single barrier in for the hole, the hole can 
tunnel through the barrier and be transmitted. If there 
are multiple barriers, there is the possibility of resonant 
tunneling through quasi-bound hole-states on the chain. 
These resonant states will appear as discrete peaks in the 
transmission probability as a function of energy. Loosely 
speaking, each quasi-bound state corresponds to a partic- 
ular single-particle cigenstate delocalized over the entire 
chain. 

For this consider a simple Hiickcl chain connected to 
left and right electrodes at the left and right most sites: 

H m = Yl + E 'XX + X) ( 2 ) 

3 Mj 

where is the hopping integral, cj and Cj are fermion 
operators which create or remove a hole from the jth site, 
and the site energy. For metal contacts, we con- 
structed two 5x5x5 simple cubic lattices each located 
at the ends of the molecular chain. This we do to insure 
that the density of states of the contacts is properly ac- 
counted for. The DNA sequence itself was connect to the 
center of the opposing faces of the left and right clusters. 
The cluster Hamiltonians were taken to be tight-binding 
with 

H m = J2 <X>< + E 4 C) (X + X) ( 3 ) 

where the hopping integral couples nearest neighbor sites. 
Finally, the clusters and DNA strands were connected via 

V Lmc = t mc (b lC ]+bl Cj ) (4) 
V Rmc = t mc {b r c\ + blc k ) (5) 

where the subscripts indicate the connected sites between 
the bridging DNA strand and the contact. 

The various 4-base-pair chains considered here are de- 
picted in Fig. The dashed lines between chains the 
indicate non-zero off-diagonal transfer integrals between 
sites on opposing chains. Each nucleotide site has associ- 
ated with it a hole-wavefunction a site energy, and 

a transfer integral t^, to the other sites in the sequence. 
Since we are dealing with a tight-binding model, connec- 
tivity between sites plays a central role in determining 
the energetics of the system. We consider the following 
general scheme: For intrachain transport, hopping occurs 
between nearest neighboring pairs. For interchain hop- 
ping, we assume that the hole on the nth unit can hop 
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Dual Single Pairing 

Contact Contact Contacts 




FIG. 2: Types of electrode-DNA-electrode contacts consid- 
ered herein. In the last case, complementary strands are con- 
nected to the opposing electrodes. 

to an unit on the other chain one site removed, nil. 
Furthermore, we assume that interchain hops can occur 
between units of the same type, 
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(0.0032eV) 
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G„ 


- G' n±1 


(0.019eV) 


C n 


- c n±1 


(0.0007eV) 


A n 




(0.035eV) 



or between G — A' and A — T" interchain neighbors 

term t\™^ 
G n ~A' n±1 (0.021eV) (7) 
A n ^T n±1 (0.016eV) 

The logic behind this scheme is that there is better tt 
overlap between interchain neighbors displaced by one 
step than those directly adjacent. This allows rela- 
tively facile hopping to occur between neighboring in- 
terchain pairs. This scheme has been verified via quan- 
tum chemical calculations by VoityukQ and the pa- 
rameters for our model we take from a recent study 
by Lakhno, Sultanov, and Pettitt 1] who considered a 
combined hopping/super-exchange model for computing 
hole-transfer rates through such model chains. The val- 
ues of the interchain hopping terms and remaining pa- 
rameters arc available as supplementary material. Since 
all of hopping terms are positive, these can be consid- 
ered as barriers for the derealization of the hole over 
connected units. 



B. Embedding system in continuum 

To connect the cluster-DNA-cluster system to a con- 
tinuum, we consider the system as if it were embedded in 



between two semi-infinite continua acting as reservoirs. 
Once we assume this, we can partition the state space 
of the total system into three domains, Qlu and Q Ra 
which span the states of the left and right electrodes and 
P which spans the electronic states of the bridging DNA 
molecule. Thus, the full Hamiltonian has the structure 

/ H L V mL \ 
H = VmL H m V mR (8) 
V V m R H R ) 

where the diagonal terms are the Hamiltonians for the 
uncoupled subsystems and V m L and V m R are the cou- 
plings between the DNA and the left and right electrodes. 
There is considerable ambiguity in how this partitioning 
can be constructed since the metal contacts and the DNA 
bridge are in intimate contact. Here, we take a somewhat 
middle ground and define our partition between system 
and reservoir by including the inner three layers of atoms 
in the contacts in the "molecular" subspace in order to 
properly include the surface density of states and the 
electronic interaction between the DNA and surface in 
our calculations. Once we have made this partitioning, 
we can use the Schwinger-Keldysh formalism [45] to de- 
rive the self-energy operators, El and £r, which embed 
the molecular sub-space into the continuum. As a re- 
sult, the Greens function for quantum particle moving 
through the molecular subspace is 

G(E) = {E-H + £ L + E^)- 1 . (9) 

Note that in general Si and £# are complex non-local 
operators which depend upon the scattering energy, E, 

v _ \ " VrnkYkm 1 firt\ 

E - e k + it] 

Here, V m k is the coupling between the chain and the con- 
tact and the sum is over the energy eigenvalues e/- of the 
contacts. [45j The %r\ in the denominator insures that the 
proper retarded scattering boundary conditions are en- 
forced. 

Written in tight-binding site-representation, they can 
be written in a local form 

S L ,„ = iA n (E) + A n (E) 

where A(E) shifts the site energy of the nth site. A n (E) 
has the effect of acting as an absorbing potential for the 
scattering wave function. This insures that the scatter- 
ing wave function obeys the proper asymptotic boundary 
conditions as we move away from the scattering region. 

We assume that the absorbing boundary conditions 
can be constructed as local self-energy terms within the 
subspace of the clusters by setting the site energy of any 
atom located beyond a given cut-off radius from the point 
of contact between the bridge and the cluster to be com- 
plex 

e c ,n = e + iT 
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with r = 0.2eV. In the examples that follow, we set the 
cut-off radius such that at least two atomic layers of the 
contact cluster were inside the cut-off radius and two lay- 
ers were outside the cut-off. This provided a reasonable 
compromise between having enough atoms in the cluster 
to represent the density of states of a simple cubic con- 
tinuum, and provide a reasonable absorbing boundary 
condition region for the scattering wave functions. 

C. Computing Transmission and Current 

To compute transmission and current though the DNA 
strand, we use the Landauer-Biittiker formula for relat- 
ing the current, I(V) to the transmission probability, 
T(E, V), which may depend upon the applied bias, V, 

I(V) = y J T(E, V)(f L (E - V/2) + fn(E + V/2))dEll) 

where fh,R are the Fermi distribution functions for the 
left and right hand electrodes. The transmission proba- 
bility is microcanonical rate for an electron injected on 
one end of the chain to leave on the other end as given 
by the Caroli equation, |H E3 

T{E) = ATr (Y L Gl{E)T R G r {E)) , (12) 

where T^r are the imaginary parts of the left and right 
self-energies describing the coupling of the molecular wire 
to the left and right hand continua. Finally, G r is 
the reduced Greens function describing hole propagation 
through the bridging region 

G r (E) = (E-H + £r\ (13) 

where H is the molecular Hamiltonian with £ = + Y,r 
being the self-energy operator of the hole on in the bridge. 

We can also consider I(V) in terms of the number of 
modes contributing to the current [ssj 

I=^-M(F + -F~) (14) 

where F are the potentials for the incoming and out- 
going holes and M is the number of modes. We can 
generalize this by letting the trace in T(E) be the sum 
over the probability that a given mode will contribute to 
the current at a given energy, 

T(E) = J2Pn(E)=Tr(P). (15) 

n 

Since the trace of a matrix is invariant to representation, 
we can construct the probability matrix. |44| 

p{e) = 4r^ /2 G t ( J B)r i? G(£;)r^ /2 , (16) 

and p n (E) are then the eigenvalues of P(E) which are 
bound by p n (E) £ [0, 1]. If we invert P, we have 

P- 1 = jTl 1/2 {E - H + Y,*)T^(E - H + £)r~ 1/2 (17) 



and the eigenvalues of P _1 are p" 1 . Consequently, we can 
compute the transmission and hence current by summing 
the eigenvalues of the P-matrix which we can construct 
directly from H. Any singularity which may be present 
in the Greens function no longer presents a problem 
when performing the energy integration in the Landauer- 
Biittiker equation. 

In practice, there is usually only one low-lying eigen- 
value of the P -1 matrix and this eigenvalue is typically 
on the order of unity. The remaining eigenvalues are 
typically on the order of 10 7 greater than the lowest 
eigenvalue. Consequently, we can use matrix diagonaliza- 
tion techniques optimized for finding isolated eigenvalues. 
Since only the lowest eigenvalue of P^ 1 contributes to to 
Tr[P] to any significant extent, once this eigenvalue has 
been determined, we do not need to compute the remain- 
ing eigenvalues. Furthermore, for determining N(E), we 
do not need to compute the eigenvectors of P _1 . Since 
computing P _1 and its lowest eigenvalue involves sim- 
ple matrix-vector manipulations, these tasks can be ef- 
ficiently performed on a parallel computer. Also, the 
method is universal and can be used to compute trans- 
mission through any model system given model Hamil- 
tonian and a self-energy representing the coupling of the 
model system to a continuum. 

One final comment regarding implementation of the 
Manthe-Miller algorithm is that as expressed P is a sin- 
gular matrix if Tl and Tr are localized to finite region 
and can not be formally inverted. To remedy this, we 
allow the absorbing potential to be very small in regions 
outside the absorbing region Tl^r = rj. Typically, we set 
7] to be between 10~ 7 to lCr 10 eV. 



D. Transmission through peptide chains 

We considered transmission through four model se- 
quences as shown in Fig. ^ Furthermore, to explore the 
effect of intra- vs. interchain transmission we considered 
three possible contact scenarios between the chain and 
the electrode surfaces. In the first case, the terminal sites 
on both chains are connected directly to the cluster sur- 
faces. In this case, hole injection and extraction occurs 
through both chains and can lead to destructive and con- 
structive interferences as in a two-slit experiment. Sec- 
ondly, we considered the case where only one chain, the 
one with the terminal G sites, is directly connected to the 
clusters. Here, the second chain acts as a hole acceptor 
and one expects to see only resonances coming from the 
intra-chain scattering pathways. Lastly, we consider the 
case where one chain is connected to each cluster and the 
chains are connected by their complementary base pairs. 
This scenario has a potential technological ramification 
in designing current probes which bind only to specific 
complementary sequences on the opposite cluster surface. 
In this case, current through the DNA must occur via in- 
terchain hops. This we expect to be quite weak, given 
that the interchain hopping matrix elements are roughly 
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an order of magnitude less than the intrachain hopping 
matrix elements. 

The corresponding transmissions through each se- 
quence is shown in Fig. [3] In the black curves plot- 
ted in Fig.|3-A through D, both terminal G-C units are 
"connected" to the left and right metal clusters through 
transfer integral terms. In this case, current can result 
via transfer through both opposing chains. The inter- 
chain hopping leads to constructive and destructive in- 
terferences between the two current pathways. In the red 
curves plotted in Fig.^l-A through D we consider the case 
in which only the terminal G sites are capable of trans- 
ferring holes between the contacts and the chain. Here, 
loosely speaking, scattering interferences can occur when 
the hole hops from the initial chain to the second and 
back to the first. Since this process is also present in 
the dual contact cases, we can distinguish purely intra- 
chain processes from multi-chain scattering processes by 
identifying peaks present in both spectra. 

In cases A and B we compare a —A — A— bridge to a 
— T — T— bridge. In both cases, there is a weak broad 
resonance at E = 0. However, the primary effect is to 
switch the order of the single resonance at E = 0.8eV 
with the cluster of three resonances at E = Q.5eV. It 
is also interesting to note that the gap between the sin- 
gle resonance and the cluster is almost unchanged by the 
switch. Comparing the dual contact transmission to the 
single contact transmission, we note that the single res- 
onance peak is present in both in the single and duel 
contact cases and that two of the three resonances in the 
cluster about E — Q.beV are present in the single con- 
tact case. For B the resonance at E ~ OAeV is present 
in both cases; however, in the single contact case, only 
one of the three corresponding resonances peaks from the 
duel contact case is present. In both A and B, the weak 
resonance at E = OeV disappears when there is only a 
single contact. 

In case C we have a — T — A— linkage and we can see 
clearly three resonance peaks, a doublet centered at 0.4 
eV and a single peak at 0.7 eV. In each of these cases, 
the corresponding eigenstate of the isolated molecular 
Hamiltonian is a delocalized interchain state. The other 
eigenstates are more or less localized to one chain or the 
other and do not extend fully across the bridge. Chang- 
ing the DNA-metal contact so that only the terminal 
G's are connected profoundly diminishes the intensity of 
each transmission resonance. Consequently, for the case 
of —T — A— linkages, interchain pathways appear to play 
an important role In case D, we transpose the linkage in 
case C to — A — T— . Here, we observe the same resonance 
peaks as in case C, except with uniformly greater inten- 
sity even when we change the nature of the DNA-metal 
contact. 

Next, we consider the sensitivity of the transmission 
to interchain fluctuations. Within our model, the in- 
terchain hopping terms represent the tunneling splitting 
between two neighboring interchain sites. Since tunnel- 
ing is highly sensitive to distance of separation, struc- 
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FIG. 3: (COLOR) Microcanonical transmission probabilities 
for hole transmission through various DNA sequences. Black 
curves correspond to transmission with both terminal G and 
C groups connected to the surface of the electrodes, Red and 
blue curves correspond to only having the terminal G group on 
one strand attached to both electrodes. For the blue curves, 
we randomly the interchain matrix elements to determine sen- 
sitivity to interchain interactions. The units for the energy 
axis are eV. 



tural fluctuations should have a profound effect on the 
the ability of a hole to tunnel between the chains. Such 
fluctuations can be included in the present model by sam- 
pling the interchain hopping matrix elements from a nor- 
mal distribution centered about some average value. For 
the blue curves, we allowed each off-diagonal interchain 
term to fluctuate about its average value by <jj = O.OleV 
and then sampled over an ensemble of 100 configurations. 
Comparing the blue and red curves in Fig|3|A., the nar- 
row resonance at 0.4eV is more or less washed out en- 
tirely by the interchain fluctuations. On the other hand, 
the two other resonances are clearly present in both the 
static case and in the fluctuating case indicating that at 
least some of the resonance features are insensitive to 
interchain fluctuations. These surviving resonances are 
certainly due to intrachain transmission. For the other 
three cases, the main resonance features are more or less 
unaffected by the interchain fluctuations. 

In the last case, we connect one strand to the left con- 
tact and the complimentary strand on the other such 
that current can only pass through the sequence through 
interchain hopping. Such a scenario could be useful in 
developing a sequence specific probe since current can 
only flow if the complementary strands match. Unfortu- 
nately, however, our calculations indicate essentially zero 
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transmission through any of the strand combinations con- 
sidered above. This is somewhat surprising considering 
that in each case at least one eigenstate of the isolated 
chain Hamiltonians arc dclocalized over the entire chain 
and have some appreciable amplitude on both of the ter- 
minal contact sites. 



III. DISCUSSION 

These results are both encouraging and discouraging 
for the potential utility of using current /voltage probes 
to assay sequence. They are encouraging since clearly 
sequence has a very profound impact on the hole trans- 
mission probability. If one assumes that a given chain 
can be decomposed into a series of A-T potential barri- 
ers separated by a series of G's and that net transmis- 
sion probability through the chain can computes as the 
product as a series of incoherent transition probabilities 
at a given energy, one could effectively engineer a DNA 
wire with specific voltage turn-on characteristics. For 
instance, having paired -A- A- or -T-T- sections as in 

5' -G-A-A-G G-T-T-G-3' 

would lead to a turn-on voltage of about 0.4 eV since 
both chain segments A and B have a series of resonances 
starting at this energy. Placing a -G-T-A-G- segment on 
the chain as in 

5' -G-A-A-G-G-T-A-G-G-T-T-G-3' 

would lead to a slightly higher turn-on voltage since the 
transmission probability for the -G-T-A-G- (sequence D) 
has a series of resonance peaks at «0.45 eV and this is in 
the middle of the cluster of resonances for sequences A 
and B. Moreover, transposing the A and T to -G-A-T-G- 
(sequence C) as in 

b'-G-A-A-G-G-A-T-G-G-T-T-G-tf 

will result in a similar turn-on voltage for the current, 
however, the current will be roughly half that of when 
we have -G-T-A-G-. 

What is also surprising is that the transmission is dom- 
inated by intrachain transport rather than a combination 
of intra- and interchain hopping. This is surprising given 
that the interchain hopping terms generally less positive 
then the intra-chain hopping matrix elements. A posi- 
tive hopping term between sites implies that the lower 
energy eigenstates will be the more localized than the 



higher energy states. The greater the hopping term, the 
more localized the lower energy states become. This is 
the opposite case from conjugated polyenes where the 
hopping term is negative leading delocalized low energy 
states and resonance stabilization. Consequently, the de- 
localized interchain states should be lower in energy than 
the delocalized intrachain states. The problem then is 
that they are insufficiently delocalized over the chain to 
provide a conduit from one end of the bridging molecule 
to the other. 

Finally, we note that the transmission only occurs 
above a certain threshold energy determined by the scat- 
tering resonances of the bridging system. Equally impor- 
tant is density of available states in the regions on ether 
side of DNA bridge. Here, we took the contacts to be 
three-dimensional simple-cubic solids with a fairly high 
density of states in the energies considered. However, in 
a real DNA molecular wire, the bridging region is more 
like 

5' - G G - (bridge) - G G - 3' 

and one should use a sequence of G's in building the 
density of states and the self-energies. If we consider a 
sequence of 5 — G — • — G— to be a quasi-one dimensional 
chain with hopping integral 7 between neighboring sites, 
then the dispersion relation for the chain is simply 

£/c = 27Cos(/c) 

with k in units of h/2a and a being the lattice spacing. 
Consequently, if the bandwidth 27 is less than the res- 
onance energy of the bridge, the transmission through 
the bridge will not occur since the density of incoming 
states at that energy is zero. On the other hand, if 7 is 
large, and one has a wide band for the holes in the gua- 
nine chains, then transmission through the bridge can 
occur. This suggests that the conductivity properties of 
a DNA chain are highly dependent upon both the base- 
pair sequence of the DNA chain and the mobility of holes 
between neighboring G sites. 
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